Rapid Adiabatic Passage from an Atomic to a Molecular Condensate 
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We examine collective magnetoassociation of a Bose- Einstein condensate (BEC), focusing on rapid 
adiabatic passage from atoms to molecules induced by a sweep of the magnetic field across a wide 
10 G) Feshbach resonance in Rb. This problem raises an interest because strong magnetoas- 
sociation is expected to favor the creation of molecular-dissociated atom pairs over the formation 
of molecular BEC [Javanainen and Mackie, Phys. Rev. Lett. 88, 090403 (2002)]. Nevertheless, the 
conversion to atom pairs is found to depend on the direction of the sweep, so that a system initially 
above threshold (open dissociation channel) may in fact give efficient conversion to molecules. 



Photoassociation occurs when a laser is on resonance 
with a transition from the coUisional state of an atom pair 
to the bound state of a molecule. The analogous process 
of magnetoassociation, which involves a Zecman-shifted 
molecular level on Feshbach resonance with a colliding 
atom pair, has an identical formalism, and intuition de- 
veloped in one case is applicable to the other. Specifically, 
the probability of photoassociation in a Bose-Einstein 
condensate (BEC) was found to approach unity [1], which 
opened the door to the possibility of using coherent mag- 
netoassociation [2] or photoassociation [3,4] of already- 
Bose-condensed atoms to create a molecular BEC. 

Early theory of coherent association was implicitly 
based on few-mode models that neglect transitions to 
noncondensate modes [2-4]. Such "rogue" [5,6], or un- 
wanted [7,8], transitions to noncondensate modes occur 
because the dissociation of a zero-momentum (k — 0) 
BEC molecule need not take the atoms back to the zero- 
momentum atomic condensate, but may just as well end 
up with two atoms with opposite momenta (±k). Since 
the coherent condensate coupling scales like the square 
root of the laser intensity (Feshbach-resonance width) 
and the noncondensate coupling scales like the intensity, 
rogue dissociation is expected to play a dominant role in 
strongly coupled atom-molecule systems [5,6]. 

Pioneering experiments [9] with photoassociation of 
*^Rb BEC were found to be just on the verge [5] of coher- 
ent atom-molecule conversion. Next-generation Na [10] 
and ''Li [11] experiments were aimed at the strongly inter- 
acting regime, and probed the potential photoassociation 
rate limit. Meanwhile, experiments on magnetoassocia- 
tion in *^Rb have led to dramatic losses of BEC atoms 
for a sweep of the magnetic field across resonance [12], 
a collapsing condensate with bursts of atoms emanating 
from a remnant condensate [13], a counterintuitive de- 
crease in condensate losses for an increasing interaction 
time [14], and collective burst-remnant oscillations [15]. 

Rogue dissociation is the key to an overall understand- 



ing of these experimental results. First, although the 
experiments remain inconclusive on this score [10,11], 
a rate limit is to be expected when the conversion to 
rogue-dissociated atom pairs dominates over the forma- 
tion of molecular condensate [6,7]. Similarly, the losses 
for across- resonance sweeps [12], and the counterintu- 
itive losses [14], can be viewed as rapid adiabatic passage 
from BEC to rogue-dissociated atoms pairs [16]. The 
remnant-burst oscillations [15] are identified as Ramsey 
fringes in the evolution between an atomic condensate 
and a molecular condensate dressed by noncondensate 
atom pairs [16-19]. 

There remains the matter of the formation of a molec- 
ular condensate. In magnetoassociation experiments, the 
sum of the remnant and burst atoms does not entirely ac- 
count for the initial condensate, and the missing atoms 
are roughly consistent with calculations [16-19] of the 
number of molecules formed, thus indirectly confirming 
the realization of molecular BEC. Either way, the (mea- 
sured) calculated fraction of (atom loss) molecular con- 
densate is small (^ 10%). 

The question is whether an improvement of the *^Rb2 
conversion efficiency is possible. Given the existing re- 
sults [3,5-8,11,13-16], the immediate answer is no. The 
155 G magnetoassociation resonance in ^^Rb is excep- 
tionally strong, and rogue dissociation is apparently 
prominent. Indeed, the qualitative agreement between 
experiment [12] and theory [16] indicates that the obvi- 
ous solution, a ~ 1 ms sweep of the magnetic field across 
resonance, merely causes rapid adiabatic passage from 
the initial atomic condensate to rogue-dissociated atoms, 
with very little molecular BEC formed. The present Let- 
ter reports our startling observation that, simply by re- 
verting the direction of the sweep of the magnetic field 
from the experiments of Ref. [12], rapid adiabatic pas- 
sage [3,5,20] may nonetheless efficiently convert from an 
atomic to a molecular condensate. 

We continue to refine our mean-field model for pho- 
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toassociation and Feshbach resonance [6,16] in an at- 
tempt to extract as much quantitative accuracy out of it 
as possible. Since Bosc enhancement favors the return of 
dissociated atoms back to the molecular BEC, we ignore 
all noncondensate modes of the molecule. The model 
thus includes an atomic condensate [described by the 
complex amplitude a], a molecular condensate [/3], and 
noncondensate atom pairs with zero total momentum 
[C(e) cx (flpfl-p), e = ?ik^/m]. As new features we now 
incorporate possible Bose enhancement in the conversion 
of molecules to noncondensate atom pairs, conversion of 
the "anomalous" atom pair amplitudes [(apa_p)] to usual 
atomic probabilities [P{e) = (a^Op)], and a gradual high- 
energy cutoff of the coupling of molecules to atom pairs 
[/(e)]. The system to be solved reads 



la = = pa , 

iP = 6oP — aa 

V 2 



(la) 



iC{e)=eC\e) 



/ de^V2P(e) + l/(e)C(e), (lb) 



The "bare" energy of the bound molecular state refer- 
enced to the dissociation threshold is hSo, the coupling 
between atomic and molecular condensates is i7, and the 
coupling for dissociation of molecules into noncondensate 
atom pairs is ^. The couplings are related in a density- 



dependent manner: ^ = Q/{2y/TTL 



The atom pair amplitude C(e) is scaled so that Eqs. (1) 
preserve the norm + + j de |C(e)|^ = 1. Also new 
is an average over the density profile of trapped atoms, 
which is identical for all three-dimensional Gaussian dis- 
tributions with the same peak density pQ. 

The energy profile /(e) is the most significant new ad- 
dition. The default shape is such that /(e) starts as 
/ = 1 at e = 0, and falls off smoothly with energy over 
a characteristic scale hAe. Physically, this high-energy 
(-momentum) cutoff accounts for the fact that the atom- 
molecule coupling is not a contact interaction, but has a 
nonzero range. We use simple Fourier transform meth- 
ods [16] to study the coupling-induced renormalization 
of the detuning, 6o S, and the binding energy Eb of 
the dressed molecules consisting of a coherent superposi- 
tion of the bare molecules [0\ and correlated atom pairs 
[C(e)]. The physical detuning varies with the magnetic 
field applied on the sample as M = {Bq — B)An, where 
Bq = 154.9 G is the position of a particular Feshbach res- 
onance in *^Rb, and A/x = 2.23 /xb [17] is the difference 
in the magnetic moments of the participant atomic and 
molecular states. 

To compare with close-coupled calculations [15,17], we 
have tried a number of different coupling functions /(e). 



both ad-hoc and forms originating from square-well and 
Lennard-Jones models. For the generic /(e), there is lit- 
tle qualitative difference between various choices. So far, 
the best match is for the ad-hoc profile |/(e)|^ = 9{eM — 
e) Ae/(Ae + e). The abrupt cutoff at cm = 20 (/us)-^ 
is for numerical purposes only, and was set high enough 
so that it has a minor effect on the results. The re- 
maining parameter ^ and Ae are chosen in such a way 
that the desired binding energy is obtained at two widely 
spaced magnetic fields. The parameter values used here 
are ^ = 19.32 (/is)-^* and Ae = 0.9537 (ms)-^ Com- 
parison with Rcf. [17] shows a maximum mismatch of 
about 9% around 160 G, which occurs because the close- 
coupling calculations were matched at the end points of 
the magnetic-field range considered (~ 157-162 G). 

In Ref. [17], an explicit sequence of two magnetic field 
pulses is specified with the implication that it was used 
to simulate a Ramsey-fringe experiment [15]. The frac- 
tion of condensate atoms, noncondensate atoms, and the 
sum thereof at the end of the pulse sequence are given 
as a function of the time that the magnetic field is held 
at a constant value between the pulses. However, the 
results in Ref. [17] are not consistent with the magnetic 
field during the holding period displayed in the pulse se- 
quence, 160 G, and we replace this by the value used in 
the experiment being simulated, 159.84 G. We have thus 
carried out the same calculations using the present meth- 
ods. In Fig. 1 we have copied the results from Ref. [17], 
and overlaid our calculations as plus signs. Though the 
comparison may look alarming, the apparent discrepancy 
originates from an about 10% difference in the frequency 
of the Ramsey fringes, exactly as expected based on the 
binding energy difference discussed above. Given that 
the two approaches are seemingly quite different, the 
match of the Ramsey fringes per se is excellent. 

Unfortunately, when we turn to quantitative compar- 
isons with experimental results, we still [16] receive a 
mixed message. For instance, let us revisit the experi- 
ment of Ref. [12]. The magnetic field is swept from 162 G 
to 132 G, so that the molecular dissociation channel is 
initially closed and opens when resonance is crossed. In 
Fig. 2, we borrow Fig. 2 from Ref. [12] showing the exper- 
imental fraction of atoms remaining in the atomic con- 
densate after the magnetic field sweep as a function of 
the inverse of the sweep rate, and overlaid our results as 
crosses. The agreement would be excellent, had we not 
used the peak density po = 1.1 x 10"*^^ cm^'^ (the same as 
in Fig. 1) whereas Ref. [12] reports po = 1-0 x 10^^ cm~^. 

Let us continue just the same to use the peak den- 
sity Po = 1.1 X lO^^cm"^, and ask what happens if the 
magnetic field is swept in the opposite direction, e.g., 
142 172 G. At first blush, not much; overlaying the 
lost-atom fraction on Fig. 2 would give agreement with 
the crosses for all but the fastest-sweep data point. 

However, a profound difference emerges when the 
progress of sweeps is followed in time. In Fig. 3 we show 
the fractions of atomic and molecular condensate, and 
noncondensate atoms, as a function of the magnetic field. 
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Pane (a) is for a 172 142 G sweep, pane (b) for the op- 
posite sweep, and both inverse sweep rates arc 20 /is/G. 
In a Cornish-type [12] sweep, condensate atoms arc un- 
ceremoniously converted to rogue atoms, and the molec- 
ular fraction remains neghgible [16]. An opposite sweep 
begins the same; but, as the system moves well past the 
resonance, nearly all noncondensate atom pairs are con- 
verted into molecules. This conclusion is independent of 
sweep rate, and survives the density averaging. For the 
slowest sweep rate, 85% of the atoms would have been 
converted to molecules and 1% to noncondensate atoms. 

To intuitively understand the double-adiabatic process 
that takes place here, consider a two-level system com- 
prised of an atomic BEG (state 1) and a molecular BEG 
dressed by noncondensate atom pairs (state 2). Near res- 
onance state 2 is primarily atom pairs, and so an adia- 
batic sweep of the detuning into the neighborhood of the 
threshold will necessarily convert the initial BEG into 
noncondensate atom pairs, regardless of direction of the 
sweep. The asymmetry arises because far below thresh- 
old state 2 is primarily molecular BEG. A sweep from 
above to below threshold will eventually convert the non- 
condensate pairs into molecular condensate. There is no 
such possibility for a sweep going the opposite way. 

More precisely, our hypothesis for a sweep of the de- 
tuning from an initially-open to a finally-closed dissocia- 
tion channel is as follows. The molecule may be strongly 
mixed with the atom pairs; but, the zero-energy atomic 
condensate holds the initial probability, and is off res- 
onance from the dressed eigenstates with any substan- 
tial molecular component. As the detuning is swept, 
eventually eigenstates with a significant molecular ampli- 
tude are sufficiently close to resonance with the atomic 
condensate for transitions to take place. Transitions 
go mainly to the dressed molecule, the negative-energy 
eigenstate of the system. Gontinuing the sweep then gives 
adiabatic conversion from dressed to bare molecules. 
When the detuning is swept in the opposite direction 
(say, 162 132 G), so that a dissociation channel opens 
for molecules after the resonance, the initial events are 
similar as in a sweep going the other way. However, af- 
ter the resonance is crossed, the bare molecules will not 
separate as an eigenstate of their own but remain diluted 
by the atom pairs. 

Whereas detuning-swccp conversion in a two-mode 
model is independent of sweep direction, strong rogue dis- 
sociation implies that transfer to molecules only happens 
for a sweep in which the molecular dissociation chan- 
nel closes after the resonance. Past rate theory reason- 
ably surmises that, after the dissociation channel closes, 
the molecules should not dissociate anymore [3]; however, 
while this serves as a convenient mnemonic, our present 
result is about coherent evolution, and at least directly is 
not about rates or a closing dissociation channel. 

While excellent quantitative agreement with the 
Ref. [17] calculations is found, we are experiencing persis- 
tent problems when trying to match experiments [12,14]. 
Missing are both quenching of the molecules and coUi- 



sional decoherence of atom pairs. No quantitative data 
on such processes are available at this time, but we have 
been experimenting by adding relaxation mechanisms 
into our calculations. With a hand-picked model for colli- 
sions, an adequate agreement with the experimental data 
in Fig. 2 is reached without adjusting the density in the 
calculations. The problem is that a universal model for 
collisions that would lead to an agreement between the- 
ory and all of the experimental data at once is elusive. 

Nevertheless, we propose making a ^^Rb2 conden- 
sate simply by asymmetric adaptation of the single- 
pulse experiments of Glaussen et al. [14]: form a BEG 
at 166 G; sweep fast (^ 1 G/fj,s) to 142 G, losing 
a negligible number of atoms (^ 1 %); sweep slowly 
(~ 5 X 10-3 G//XS) back to 172 G, efficiently con- 
verting the atoms to molecules. For modest densities 
((p) ~ IQ-^^cm-^), ro-vibrational quenching at a rate 
comparable to Na [2] should be tolerable on millisecond 
timescales [21], and the efficient production of a molecu- 
lar BEG ought to be within easy reach of experiments. 

We have come full circle. In our earliest publication on 
coherent photoassociation (and, by formal equivalence, 
on magnetoassociation) of a BEG, we predicted rapid adi- 
abatic passage from an atomic to a molecular condensate 
induced by a sweep of the detuning S across the resonance 
[3]. Those arguments were based on a two-mode model, 
i.e., atomic and molecular condensates only. In the in- 
terim it has become evident [5.6] that rogue dissociation 
to noncondensate atoms is a dominant factor in the dy- 
namics whenever fl ~ LOp, a, condition met by two orders 
of magnitude in the JILA experiments [12-15]. Hence, 
rapid adiabatic passage was not a leading possibility for 
creating a ^^Rb2 molecular condensate. Now it has made 
a comeback with a vengeance. 

While preparing this manuscript, we became aware of 
similar Na work by Yurovsky and Ben-Reuven [22] . 
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FIG. 2. Experimental (•, error bars) [12] and theoretical 
(x) fraction of atoms remaining after the magnetic field is 
swept linearly from 132 G to 162 G, versus the inverse sweep 
rate. The theoretical (experimental) peak density of the con- 
densate is po = 1.1 X lO" cm"^ (po = 1-0 x 10^^ cm"^). 
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FIG. 1. Ramsey fringes in atom-molecule conversion, as 
calculated herein {+) and in Ref. [17] (o,-, — ,□). From the 
bottom up, we have the number of noncondensate and con- 
densate atoms, and the sum thereof. The Ramsey pulses 
are as in Fig. 2 of Ref. [17], except that the hold field was 
159.84 G. The peak density is po = 1-1 x 10^^ cm"^. 
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FIG. 3. Fraction of atoms in atomic condensate (dotted 
line), in noncondensate atoms (dashed line), and in the molec- 
ular condensate (solid line) present while the magnetic field 
is swept by 1 G per 20 ps in the direction indicated by the 
arrow in each pane. Molecules are too few to resolve in the 
upper pane. There is no density average; the fixed density 
p = 3.9 X 10^^ cm~^ is what the average density would be for 
the peak density po = 1.1 x 10^'^ cm~'^. 
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